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The bound state energy eigenvalues for the two-dimensional Kepler problem are found to be 
degenerate. This "accidental" degeneracy is due to the existence of a two-dimensional analogue of 
the quantum-mechanical Runge-Lenz vector. Reformulating the problem in momentum space leads 
to an integral form of the Schrodinger equation. This equation is solved by projecting the two- 
dimensional momentum space onto the surface of a three-dimensional sphere. The eigenfunctions 
are then expanded in terms of spherical harmonics, and this leads to an integral relation in terms 
of special functions which has not previously been tabulated. The dynamical symmetry of the 
problem is also considered, and it is shown that the two components of the Runge-Lenz vector in 

C^l ' real space correspond to the generators of infinitesimal rotations about the respective coordinate 

{^J axes in momentum space. 
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I. INTRODUCTION 



A semiconductor quantum well under illumination is a quasi-two-dimensional system, in which photoexcited elec- 
trons and holes are essentially confined to a plane. The mutual Coulomb interaction leads to electron-hole bound 
states known as excitons, which are extremely important for the optical properties of the quantum well. The rel- 
■ ative in-plane motion of the electron and hole can be described by a two-dimensional Schrodinger equation for a 
single particle with a reduced mass. This is a physical realization of the two-dimensional hydrogenic problem, which 
I/"} ■ originated as a purely theoretical construction JD. An important similarity with the three-dimensional hydrogen 
atom is the "accidental" degeneracy of the bound state energy levels. This degeneracy is due to the existence of the 
, quantum- mechanical Runge-Lenz vector, first introduced by Pauli in three dimensions, and indicates the presence 
' of a dynamical symmetry of the system. 

The most important study relating to the hidden symmetry of the hydrogen atom was that by Fock in 1935 
||. He considered the Schrodinger equation in momentum space, which led to an integral equation. Considering 
negative energy (bound-state) solutions, he projected the three-dimensional momentum space onto the surface of a 
four-dimensional hypersphere. After a suitable transformation of the wavefunction, the resulting integral equation 
was seen to be invariant under rotations in four-dimensional momentum space. Fock deduced that the dynamical 
symmetry of the hydrogen atom is described by the four-dimensional rotation group SO (4), which contains the 
geometrical symmetry SO (3) as a subgroup. He related this hidden symmetry to the observed degeneracy of the 
. £h , energy eigenvalues. 

Shortly afterwards, Bargmann [Q made the connection between Pauli's quantum mechanical Runge-Lenz vector 
and Fock's discovery of invariance under rotations in four-dimensional momentum space. Fock's method was also 
extended by Alliluev Q to the case of d dimensions (d ^ 2) . A comprehensive review concerning the symmetry of the 
hydrogen atom was later given by Bander and Itzykson |j, including a detailed group theoretical treatment and 
extension to scattering states. 

Improvements in semiconductor growth techniques over the subsequent decades, which enabled the manufacture 
of effectively two-dimensional structures, led to a resurgence of interest in the two-dimensional hydrogen atom. The 
Runge-Lenz vector for this case was defined for the first time S, and real-space solutions of the Schrodinger equation 
were applied to problems of atomic physics in two-dimensions iB . 

Recent studies have focused on diverse aspects of the hydrogenic problem. The d-dimensional case has been 
reconsidered, leading to a generalized Runge-Lenz vector (see 10 and references therein). The algebraic basis of the 
dynamical symmetry has also been given a thorough mathematical treatment JH], [l^] . 

In the present work we return to the two-dimensional problem, and use the method of Fock to obtain a new 
integral relation in terms of special functions. The dynamical symmetry of the system is also considered, and a new 
interpretation of the two-dimensional Runge-Lenz vector is presented. 
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II. PROBLEM FORMULATION 
A. Preliminaries 



The relative in-plane motion of an electron and hole, with effective masses m e and m/,, respectively, may be treated 
as that of a single particle with reduced mass p = m e mh/(rn e + mh) and energy E, moving in a Coulomb potential 
V{p). The wavefunction of the particle satisfies the stationary Schrodinger equation 
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*(p) - E9(p), (1) 



where (p,4>) are plane polar coordinates. Note that excitonic Rydberg units are used throughout this paper, which 
leads to a potential of the form V(p) = —2/ p. 

The eigenf unctions of Eq. (|l|) are derived in Appendix [a|. It is well known that the bound state energy levels are 
of the form 

E = -—±— S) n = 0,l,2,... (2) 
(n + 1/2) 

where n is the principal quantum number. Notably, Eq. (^) does not contain explicitly the azimuthal quantum number 
m, which enters the radial equation (see Appendix Eq. (|A4|) ). Each energy level is (2n + l)-fold degenerate, the 
so-called accidental degeneracy. 

It is convenient to introduce a vector operator corresponding to the z-projection of the angular momentum, L z = 
e z L z , where e z is a unit vector normal to the plane of motion of the electron and hole. We now introduce the 
two-dimensional analogue of the quantum-mechanical Runge-Lenz vector as the dimensionless operator 

2 

A = (q x L z - L z x q) - -p, (3) 

P 

where q = — iV is the momentum operator. Note that A lies in the plane and has Cartesian components A x and A y . 
L z , A x and A y represent conserved quantities and therefore commute with the Hamiltonian: 

[H,L z ] = [H,A x ] = [H,Ay}=0. (4) 

They also satisfy the following commutation relations: 

[L z ,A x ] = iA y , (5) 
[L z , A y ] = -iA x , (6) 
[A x ,A y ] = -4iL z H. (7) 

B. Derivation of energy eigenvalues from A 

The existence of the non-commuting operators A x and A y , representing conserved physical quantities, implies that 
the Runge-Lenz vector is related to the accidental degeneracy of the energy levels in two dimensions jl3) . We now 
present a simple interpretation of the hidden symmetry underlying this degeneracy. 

For eigenfunctions of the Hamiltonian we can replace H by the energy E, and defining 

A' = — (8) 

2y^E 



we obtain the new commutation relations: 



[L z ,A' x ]=iA! y , (9) 
[L z , A' y ] = -iA' x , (io) 
{A' x ,A' v ] = iL z . (ii) 
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If we now construct a three-dimensional vector operator 

J = A' + L Z , (12) 

then the components of J satisfy the commutation rules of ordinary angular momentum: 

[Jj, J k ] = itjklJh (13) 

where ejki is the Levi-Civita symbol. 

Noting that A' • L z = L 2 ■ A' = 0, we have 

J 2 = (A' + L 2 ) 2 = A' 2 +f4, (14) 

where the operator J 2 has eigenvalues j(j + 1) and commutes with the Hamiltonian. 

We now make use of a special expression relating A 2 and L 2 , the derivation of which is given in Appendix 

A 2 = (4L 2 + 1) + 4. (15) 

Substituting in Eq. ( [l4|) and again replacing H with E, we obtain 



J 2 = -^ 



£(4L 2 + l)+4 +L 2 . (16) 



Because [H, J 2 ] = 0, an eigenfunction of the Hamiltonian will also be an eigenfunction of J 2 . Operating with both 
sides of Eq. (|l^) on an eigenfunction of the Hamiltonian, we obtain for the eigenvalues of J 2 : 

jC7 + D = -(i + |)- ( 17 ) 

Rearranging, and identifying j with the principal quantum number n, we obtain the correct expression for the energy 
eigenvalues: 

E = ~ ( * , 2 , n = 0,l,2,... (18) 
(n + 1/2) 

Note that the z-component of J is simply L z . Recalling that the eigenvalues of L z are denoted by m, there are 
(2j + 1) values of m for a given j. However, as j — n, we see that there are (2n + 1) values of m for a given energy, 
which corresponds to the observed (2n + l)-fold degeneracy. 

III. FOCK'S METHOD IN TWO DIMENSIONS 

A. Stereographic projection 

The method of Fock ||, in which a three-dimensional momentum space is projected onto the surface of a four- 
dimensional hypersphere, may be applied to our two-dimensional problem. We begin by defining a pair of two- 
dimensional Fourier transforms between real space and momentum space: 

$(q) = J ^{ P y*Pdp, (19) 
*(P) = p^a / $(q)e^ q p dq. (20) 

We shall restrict our interest to bound states, and hence the energy E = will be negative. 
Substitution of Eq. (|2(]) in Eq. (|l|) yields the following integral equation for ^(q): 
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The two-dimensional momentum space is now projected onto the surface of a three-dimensional unit sphere centered 
at the origin, and so it is natural to scale the in-plane momentum by qo. Each point on a unit sphere is completely 
defined by two polar angles, 9 and <f>, and the Cartesian coordinates of a point on the unit sphere are given by 

Zqoqx 
q 2 + q^ 

q 2 + qo' 



u x = sin cos <p = — = , (22) 

T + qo 2 

u y = sin sin 4> = ^ qo<ly - , (23) 



q 2 - qo 2 



q +qo 

An element of surface area on the unit sphere is given by 



Uz =cos0= I " u 2 . (24) 



<ft> : smH<1H,k, = ( o 2g0 J ] dq, (25) 



and the distance between two points transforms as: 

2go 

(<7 2 + 9 )' /2 (^ 
If the wavefunction on the sphere is expressed as 



u-u'| = ^ TTolq-q'l- (26) 



x(u) = ^(^M *(q), (27) 



(28) 



./3o V 2 «o 

then Eq. ( pl| ) reduces to the simple form: 

XW 2-KqJ |u-u'| 

B. Expansion in spherical harmonics 

Any function on a sphere can be expressed in terms of spherical harmonics, so for x( u ) we have 

oo / 

X (u)=^^i lm r(^), (29) 

1=0 m=-l 

where Y l m (9, cf>) are basically defined as in fl4|| : 

where P„ '(coaff) is an associated Legendre function as defined in [|l5). The constant c\ m is an arbitrary "phase 
factor". As long as |c/ m | 2 = 1 we are free to choose q m , and for reasons which will become clear we set 

Clm = H) H - (31) 

The kernel of the integral in Eq. ( |28| ) can also be expanded in this basis as |lq| : 

i oo A . 

1 1 A=0p=-A 

Substituting Eqs. @ and (f|) into Eq. (|f) we have 

OO 2 



E E A >™Yi r 

— EE E E \i^—A hmi Y^{e\$)Y™\6A)YC*^^) d ^- ( 33 ) 



Z^O m=— i 

oo oo Zi i 



^1—0 /2— mi — — Hi m,2 ——l2 
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We now make use of the orthogonality property of spherical harmonics to reduce Eq. (|3^) to the following: 

00/ ~ 00 1 1 1 

E E 4mir(M) = -E E 2i7TT Aimi ^ 1( ^ 0) - (34) 

i— m— — / Omi — — ?i 

Multiplying both sides of Eq. ( |34| ) by Y^™ 0) and integrating over dfl gives 

2 

A-nm' — 77; ; TZ-^nm'i (35) 

g (2n + lj 

where we have again used the orthogonality relation for spherical harmonics. The final step is to rearrange for qq 
and identify the index n with the principal quantum number. This enables us to find an expression for the energy in 
excitonic Rydbergs: 

E=-ql = — 1 2 , n = 0,1,2,... (36) 
[n + 1/2) 

This is seen to be identical to Eq. (||). 

For a particular value of n, the general solution of Eq. (p8|) can be expressed as 

n 

X „(u)= E A nm Y™{e,4>). (37) 

m— — n 

Each of the functions entering the sum in Eq. ( p7| ) satisfies Eq. (p8|) separately. So, for each value of n we have (2n+ 1) 
linearly independent solutions, and this explains the observed (2n + l)-fold degeneracy. 

We are free to choose any linear combination of spherical harmonics for our eigenfunctions, but for convenience we 
simply choose 

Xnm(u) = A nm Y™(9,q>). (38) 
If we also require our eigenfunctions to be normalized as follows: 

J | x (u)| 2 ^=^ J £-^|$(q)| 2 dq = J \*(p)\* d p = l i (39) 

then Eq. ( f38| ) reduces to 

X „ m (u) = 27rY n ro (0,<£). (40) 



Applying the transformation in Eq. fl27j), we can obtain an explicit expression for the orthonormal eigenfunctions 
of Eq. fel 



*„ m (q) = Cnmxh ^J™ 1 ^ ( ^ 4 H («» 0)e im *, (4) i 

(n+ |m|)! \9 +9o/ 



where we have used the fact that q = (n + 1/2) 1 , and 9 and (j> are defined by Eqs. (p3)-(|2 



C. New integral relations 

To obtain the real-space eigenfunctions ^f(p) we make an inverse Fourier transform: 

1 f 1 f 2n f°° , 



(42) 



where 4>' is the azimuthal angle between the vectors p and q. However, if we now substitute Eq. (|4l]) into this 
expression we have to be careful with our notation. The angle labeled (p in Eq. (fll]) is actually related to 4>' via 

^ = + tf> P , (43) 
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where 4> p is the azimuthal angle of the vector p, which can be treated as constant for the purposes of our integration. 
Taking this into account, the substitution of Eq. ( |4l| ) into Eq. ( (42|) yields: 



Vu>)= °" m J— l TO ^ ! e im <t>P f°° ( 2q ° \ 1 p\™\( c0s n) e *(™<t>'~ip™s4>') adadd' (44) 

* (P) (2tt)3/2 V(n+H)! A io U 2 + ?oV " ^ j 70 ' 



From Eq. we obtain 



pM(coB0)=J»M(*^j | (45) 



and we use the following form of Bessel's integral p 



2tt 

e W-iP"»*') dcj)' = 2ir(-i) m J m (qp), (46) 



where J m {<ip) is a Bessel function of the first kind of order m. Substituting Eqs. (45) and (]46| ) into Eq. ( |44| ) leads to 

"f(ra)^(?Tp)^"* (47) 



'(n- 


\m\)\ 


(n + 


\m\)\ 



We now make a change of variables, x = qop and y = (? 2 /<? 2 , so that Eq. ( |47| ) becomes 
where we have used the fact that |16|] : 

p|m 



(ftt) = (^f) ■ (49) 

If we now equate the expression for in Eq. ( f48| ) with that derived in Appendix [A|, we obtain the following: 

/ i\r.+m/ -Nm [°°p\m\( 1 -y \ J m{Xy/y) (2x)\ m \e- X 2\ra\ , 



The value of c nm chosen earlier in Eq. ( |31| ) ensures that both sides of Eq. (|5C|) are numerically equal. If we restrict 
our interest to m ^ then the relation simplifies to 

1-lA Jmixy/y) _ (-l)"(2x) m e -% 2m 



P n ( TTf ) (TT#^ * = n+i/2 ^"^ n, m = 0, 1, 2, . . . ; m < n. (51) 



As far as we can ascertain, this integral relation between special functions has not previously been tabulated. For 
n, m = we recover the known integral relation Ea |: 

IV. DYNAMICAL SYMMETRY 
A. Infinitesimal generators 



Consider now a vector u from the origin to a point on the three-dimensional unit sphere defined in Sec. [II A. If 
this vector is rotated through an infinitesimal angle a in the (u x u z ) plane, we have a new vector 



u' = u + Su, 



(53) 
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where the components of u are given in Eqs. (E2)-(EJ), and 



Su = a e y x u. 



(54) 



This rotation on the sphere corresponds to a change in the two-dimensional momentum from q to q'. The Cartesian 
components of Eq. (p3) are then found to be 



2qoq' x 2q q x 



q 12 + qo r + qt> 

Zqoq'y 2q Q q y 

q' 2 + q 2 ~ q 2 + q 2 ' 

/2 2 2 2 

q -% q - % 



g 2 g 2 o 

\ 2 + qV 



2qoq x 



2 q' 2 + q 2 q 2 + ql q 2 + qV 

where q 2 = q 2 + q 2 . 

After some manipulation we can also find the components of <5q = q' q: 

q 2 - ql - iql 



Sq x = a- 
5q y = —a- 

The corresponding change in $(q) is given by 

'q 2 -ql-^q 2 x 9 q x q y d 



2qo 
q_xqy 

qo 



<5$(q) - 



(q 2 + qlf' 2 



We can write this as 



2<Zo 



<J*(q) 



where the infinitesimal generator is given by 

i 



Ax — 



{q 2 +qlY' 2 



dq x qo 



-— aA x $(q) 

^qo 



d 



(q'+q 2 ) 3 / 2 ^) 



d 



(q z - q 2 - 2qi)— - 2q x q y — + q$) 



dq x 



dq y 



-,2\3/2 



We now make use of the following operator expression in the momentum representation: 

p = e x x + e y y = iV q , 

and the commutation relation 

[p, /(q)]=iV g /, 

to derive a more compact expression for A x : 

A x = (q 2 - ql)x - 2q x {(\ ■ p) - 3iq x . 
By considering an infinitesimal rotation in the (u y u z ) plane we can obtain a similar expression for A y : 

A = (q 2 - ql)v - 2 %(<i • p) - aq y - 



(55) 
(56) 
(57) 



(58) 
(59) 

(60) 
(61) 

(62) 

(63) 
(64) 

(65) 
(66) 



These expressions operate on a particular energy eigenfunction with eigenvalue — q^. If we move the constant — q^ 
to the right and replace it with the Hamiltonian in momentum space, Ti: 

A x = q 2 x + xTL - 2g x (q • p) - 3iq x , (67) 
A y = q 2 y + yH-2q y {( l -p)-3iq y , (68) 

then A x and A y can operate on any linear combination of eigenfunctions. 



8 



B. Relation to Runge-Lenz vector 

Recall the definition of the two-dimensional Runge-Lenz vector in real space: 

2 

A = (q x L 2 - L z x q) p. (69) 

P 

Using 1L Z = p x q, and the following identity for the triple product of three vectors: 

ax (b x c) = (a-c)b- (a-b)c, (70) 
we can apply the commutation relation [ p, q] = i to rewrite Eq. ( |69|) in the form: 

A = q 2 p + p(V-^ - 2q(q • p) - 3zq. (71) 

If we now return to the expression for the real-space Hamiltonian in Eq. ([!]) , it is apparent that we may substitute 

q 2 — 2/p = H (72) 

in Eq. @ to yield 

A = q 2 p + pH — 2q(q • p) - 3iq. (73) 

Comparing this with Eqs. ( |67| ) and (|68|), it is evident that the two components of the Runge-Lenz vector in real space 
correspond to the generators of infinitesimal rotations in the (u x u z ) and (u y u z ) planes. 

V. CONCLUSION 

We have shown that the accidental degeneracy in the energy eigenvalues of the two-dimensional Kepler problem 
may be explained by the existence of a planar analogue of the familiar three-dimensional Runge-Lenz vector. By 
moving into momentum space and making a stereographic projection onto a three-dimensional sphere, a new integral 
relation in terms of special functions has been obtained, which to our knowledge has not previously been tabulated. 
We have also demonstrated explicitly that the components of the two-dimensional Runge-Lenz vector in real space 
are intimately related to infinitesimal rotations in three-dimensional momentum space. 

APPENDIX A: SOLUTION OF REAL-SPACE SCHRODINGER EQUATION 

We apply the method of separation of variables to Eq. (0), making the substitution 

M>(p) = R(pM^). (Al) 
Introducing a separation constant m 2 , we can obtain the angular equation 

d 2 $ 



, , m 2 $ = 0, (A2) 



$(0) = — Le™*. (A3) 



with the solution 

<S>(6) = 

/2tc 

The corresponding radial equation (with E = —q^) is 

d 2 R IdR (1 , m 2 \ , k . 

77 + -^+ [--11 - — )R=0. A4 
ap z p dp \p p z J 

We make the substitution 

R(p) = Cp^e- qop w(p), (A5) 
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P— + ( 2 H + 1 - 2 <Zop) — + (2 - 2\m\q Q - q )w = 0. (A6) 



where C is a normalization constant. This leads to the equation 

d 2 w „ „ , dw 

,_ + ( 2 H +1-2^)^ 

Making a final change of variables j3 = 2qop, we obtain 

dw 



z 3 ^- + ( 2 H + 1 " + ( — - H - o ) w = °- ( A? ) 



This is the confluent hypergeometric equation |L5), which has two linearly independent solutions. If we choose the 
solution which is r egu lar at the origin, then this becomes a polynomial of finite degree if qo — (n + 1/2) _1 with 



n = 0, 1, 2, . . . Eq. (|A7| ) then becomes the associated Laguerre equation |lq| . the solutions of which are the associated 

-2|ro| to\ _ r 2|m| 



Laguerre polynomials: 

» = iX|C') = i;,'H,,,'-w). ian. 



We can now write the real-space wavefunction in the form 



*nm(p) = ^ ml e^L^ m l ml (2q oP )e" n ^, (A9) 



where the reason for the subscript on <f> is explained in Sec. Ill C| 



To normalize this wavefunction we need to make use of the integral |16| 



oo 

e 



o 



2 -P(2 g op) 2 ' m l +1 i 2 Jl|L | (2<zop)^ l :;.L|(2*p) d(2q oP ) = |^|(2n + 1). (A10) 



The normalized wavefunctions are therefore: 



= fit + H)! ^ Me ~™ L n-U^P^ P > ( AU ) 

satisfying the following orthogonality condition: 

/ *n imi (p)*n 2 m 2 (p) dp = S (A12) 



APPENDIX B: DERIVATION OF EQ. 15 



From Eq. (|3|) we have 

f 2 I 2 

A 2 = (q x L z - L 2 x q) p 

P . 

2 - 2 

= [2(q x L z ) - iq] 2 - -p • [2(q x L z ) - iq] - -[2(q x L z ) - iq] • p + 4. 
P P 

We further expand as follows: 

[2(q x L z ) - iq] 2 = 4(q x L z ) 2 - 2iq ■ (q x L z ) - 2i(q x L z ) • q - q 2 
= 4q 2 L 2 + 2q 2 q 2 = q 2 (4L 2 + 1), 



and 



-p ■ [2(q x L.) - iq] - -[2(q x L z ) - iq] • p = --(4L 2 + 1). 
P P P 



Substituting Eqs. (B2) and (B3) into Eq. (Bl) gives 



A 2 = q 2 (4L 2 + l)-±(4L 2 + l) + 4, 



(Bl) 



(B2) 



(B3) 



(B4) 
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which, from Eq. (f72f), is just 

A 2 = H(4L 2 + l)+4. (B5) 
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